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transition of the theory at a critical number of fermion flavors Nf c . Our preliminary results indicate 
that Nf c « 1.5. 
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1. Introduction 



The study of QFTs in which the ground state shows a sensitivity to the number of fermion fla- 
vors Nf is intrinsically interesting. According to approximate solutions of continuum Schwinger- 
Dyson equations (SDEs), QED3 displays this phenomenon. It is believed to be confining and 
exhibit features such as dynamical mass generation when the number of fermion flavors Nf is 
smaller than a critical value Nf c [|], §, §, 0]. Apparently, for Nf > Nf c , the attactive in- 

teraction between a fermion and an antifermion due to photon exchange is overwhelmed by the 
fermion screening of the theory's electric charge. Furthermore, over the past few years QED3 has 
attracted a lot of attention, because of potential applications to models of high T c superconductivity 
[§, |9[ |6[ [T(], [II]] It is also an interesting and challenging model field theory which is being seen as 
an ideal laboratory to study more complicated gauge field theories. 

Initial studies based on SDEs using the photon propagator derived from the leading order 1 /Nf 
expansion, where Nf is the number of fermion flavors, suggested that for Nf less than Nf c ~ 3.2 
chiral symmetry is broken [jlp. The model in the limit Nf — > Nf c is supposed to undergo an infinite- 
order phase transition [g]. Other studies taking non-trivial vertex corrections into account predicted 
chiral symmetry breaking for arbitrary Nf [§]. Studies which treat the vertex consistently in both 
numerator and denominator of the SDEs have found Nf c < °°, with a value either in agreement with 
the original study [Q], or slightly higher Nf c ~ 4.3 An argument based on a thermodynamic 
inequality predicted Nf c < § [12], a result that was later on challenged in [11]. A gauge invariant 
determination of Nf c based on the divergence of the chiral susceptibility gives Nf c ~ 2.16 [01. 
Recent progress in the direction of gauge covariant solutions for the propagators of QED3 showed 
that in the Landau gauge a chiral phase transition exists at Nf c ~ 4 |JJ|. Furthermore, a perturbative 
analysis of RG flow equations in the large-A^ limit predicts Nf C ~ 6 Jl3[|. 

There have also been numerical attempts to resolve the issue via lattice simulations of QED3. 



Once again, opinions have divided on whether Nf C is finite [ |14[ , |15| ], or whether chiral symmetry 
is broken for all Nf [[If]]. A numerical study of the quenched (Nf = 0) case has shown that chiral 
symmetry is broken fll7|]. More recent numerical results showed that chiral symmetry is also broken 
for Nf = 1, whereas Nf = 2 appeared chirally symmetric with an upper bound of 10~ 4 on the 
dimensionless condensate [18]. The principal obstruction to a definitive answer has been large 
finite volume effects resulting from the presence of a massless photon in the spectrum, which 
prevent a reliable extrapolation to the thermodynamic limit. Recent lattice simulations of the three- 
dimensional Thirring model, which may have the same universal properties as QED3, predicted 
Nf c = 6.6(1) [|l9|. In this paper we present a study of the QED3 Equation of State based on 
preliminary results extracted from lattice simulations on large lattices (up to 8O 3 ). 



2. The Model 



We are considering the four-component formulation of QED3 where the Dirac algebra is 
represented by the 4x4 matrices 70, 71 and y 2 - This formulation preserves parity and gives 
each spinor a global U(2) symmetry generated by 1, 73,75 and £7375; the full symmetry is then 
U (2Nf). If the fermions acquire dynamical mass the U (2Nf) symmetry is broken spontaneously to 
U (Nf) x U (Nf) and 2Nj Goldstone bosons appear in the particle spectrum. 
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The action of the lattice model we study is 

M(x,x) = m8 x>x ! + ^£t? j u(x)[4' jc+pVxii- jc-(i U x-[i ,J- 

This describes interactions between TV flavors of Grassmann- valued staggered fermion fields 
defined on the sites x of a three-dimensional cubic lattice, and real photon fields d x ^ defined on the 



link between nearest neighbour sites x, x+fl. Since 2 is unbounded from above, eq.(^j]) defines a 
non-compact formulation of QED; note however that to ensure local gauge invariance the fermion- 
photon interaction is encoded via the compact connection C/ ZjU = exp(i6 x ^), with U x+ jx- i x = U*„. 
In the fermion kinetic matrix M the Kawamoto-Smit phases 



Tj M (x) = (-ir^-^-' (2.2) 

are designed to ensure relativistic covariance in the continuum limit, and m is the bare fermion 
mass. 



If the physical lattice spacing is denoted a, then in the continuum limit ad — > 0, eq.( |2.1[ ) can 
be shown to be equivalent up to terms of 0(a 2 ) to 

N f 1 

S = £ W J tin ih + igAp )+m]Y> + jFuvFuv (2.3) 

7=1 4 

ie. to continuum QED in 2+1 euclidean dimensions, with \\f, y describing Nf flavors of four- 
component Dirac spinor acted on by 4x4 matrices y^, and Nf = 2N. The continuum photon field 
is related to the lattice field via VjU = agA^(x), with dimensional coupling strength g given by 
g 2 = (aj8) _1 , and the field strength F^ v = d^A v — d^A v . The continuum limit is thus taken when 
the dimensionless inverse coupling j8 — > °°. 



As reviewed in [[180, for a > in the chiral limit the lattice action ( |2.1| ) retains only a remnant 



of the \J(2Nf) global symmetry of ( |2.3| ) under global chiral/flavor rotations, namely a U(/V) (g)U(/V) 
symmetry which is broken to U(/V) either explicitly by the bare mass m ^ 0, or spontaneously by 
a chiral condensate ^ 0, in which case the spectrum contains N 2 exact Goldstone modes. 
It is expected that the symmetry breaking pattern \](2Nf) — > U(7Y/) (g> U(/Y/) is restored in the 
continuum limit, implying the existence of an additional 7./V 2 approximate Goldstone modes whose 
masses vanish as /3 — > °°. 



3. Numerical Simulations 



In this section we present results from numerical simulations performed using the standard 
Hybrid Molecular Dynamics (HMD) R-algorithm. In order to ensure that the 0{N 2 dt 2 ) system- 
atic errors {dt is the time step of the HMD trajectory) are negligible, we performed two sets of 
simulations one with dt = 0.010 and another with dt = 0.005 on 80 3 lattices with m = 0.005 and 
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Figure 1: Chiral condensate versus Nf from simulations with dt = 0.010 and dt = 0.005; j3 = 0.90, m 
0.005, and L = 80. 



j8 = 0.90. As shown in Fig. [I] the values of the chiral condensate from the two sets of simula- 
tions with Nf = 0.5,..., 3.0 agree within statistical error, implying that the systematic effects are 
negligible. 

Next we present results from simulations with Nf = 1.5 and Nf = 2.0 (Figs. ^| and [|). More 
specifically, we study the behavior of the dimensionless chiral condensate /3 2 (i//i//) versus the di- 
mensionless fermion bare mass fim on different lattice sizes. In the Nf = 1.5 case the finite size 
discrepancy between the L = 54 and L = 80 results is ~ 7%, implying that the finite size effects for 
L = 80 are small. In the Nf = 2.0 case the discrepancy is ~ 1%, implying that in this case the values 
of the condensate extracted from the 80 3 simulations should be even closer to the thermodynamic 
limit values. 

For N f = 2, we fitted the 54 3 data to 

P 2 (YY) =a Q + ai ■(pm)+a 2 -(pm) 2 . (3.1) 

The extracted value of cio = — 8 x 10~ 8 with a statistical error 5 x 10~ 7 and fit quality % 2 /DOF = 
0.9 is consistent with zero with relatively high accuracy. For Nf = 1.5, we fitted the 80 3 data to 
( fQl ) and got a = -2 x 10~ 6 with a statistical error 10~ 6 and fit quality % 2 /DOF = 4.1. This 
result is also very close to zero although the fit quality is lower than in the Nf = 2.0 case. As we 
will se below the low fit quality could be attributed to the fact that Nf = 1.5 may be close to the 
chiral phase transition. 

We also checked the effects of lattice discretization on the values of the chiral condensate by 
comparing data extracted from simulations with /3 = 0.90 and j8 = 1.20 for Nf = 0.5, ...,2.0. We 
achieved this by fixing for the two sets of simulations the physical volume (L/j3) 3 and the physical 
mass /3m. The results presented in Fig.|] show that the lattice discretization effects at j8 = 0.90 are 
small for Nf > 0.5, whereas for Nf = 0.5 there is an ~ 8% discrepancy between the values of the 
condensate at j8 = 0.90 and j8 = 1.20. 
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Figure 2: Dimensionless condensate versus dimensionless bare mass for Nf = 1.5. 
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Figure 3: Dimensionless condensate versus dimensionless bare mass for Nf = 2.0. 



Next, we fitted the data for the dimensionless condensate at different values of Nf and m, and 
lattice sizes 54 3 and 80 3 to a renormalization group inspired equation of state that includes a finite 



size scaling term [20] 



m=A({P-f5 c )+CL^)(l5 2 (Y V )y + B(P 2 (W)) t 



(3.2) 



where p = 8 - l/j8 m . The results extracted from this fit are: A = 0.0477(38), B = 0.79(2), C = 
10.7(8), Nf c = 1.52(6), 8 = 1.177(7),/? = 0.73(2). The data and the fitting functions are shown in 
Fig. Bl. These results are consistent with a relatively smooth second order phase transition. 



4. Summary 

We studied numerically the equation of state of non-compact QED3. The extrapolations to the 
chiral limit on lattices with small finite size effects show with high accuracy that the Nf = 2.0 theory 
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Figure 4: Dimensionless condensate versus Nf at 3 — 0.90 and 1 .20 with constant L/3 and 3m. 
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Figure 5: Fits of j3 2 ( v F v F) vs. Nf to a finite volume scaling form of the Equation of State. 

is chirally symmetric with an accuracy of 0(1O~ 7 ). The preliminary results extracted from fits to 
a finite volume Equation of State are consistent with a second order phase transition scenario and 
Nfc ~ 1.5. However, in order to reach a decisive conclusion on the value of Nf c and the properties 
of the chiral phase transition we have to extend our simulations to larger lattices, generate better 
statistics, and include more data points close to the continuum limit. 
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